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^^ • We study dark matter(DM) in the model with one keV scale right-handed neutrino 



Vfii and two GeV scale right-handed neutrinos i'r2,3, the z^SM. We find that one of 
the GeV scale right-handed neutrinos can have much longer lifetime than the other 



^1^' when two GeV scale right-handed neutrinos are degenerate. We show that mass and 

i-G . mixing of light neutrinos can be explained in this case. Significant entropy release 

^ I can be generated in a reheating produced by the decay of one of the GeV scale ur. 

The density of i^ri DM can be diluted by two orders of magnitude and the mixing of 

in 

en ■ i^Ri with active neutrinos is allowed to be much larger, reaching the bound from X- 

cn _ 

ly-^ , ray observation. This mixing can lead to sizeable rate of z^/ji capture by radioactive 

O 

(^ [ nuclei. The I'm capture events are mono-energetic electrons with keV scale energy 

away from the beta decay spectrum. This is a new way to detect DM in the universe. 
U ■ PACS numbers: 14.60.Pq, 13.15. +g 



Introduction: Among many models of DM candidate a seesaw type model, the i^SM, 
is of particular interests l\. In this model there is one keV scale right-handed neutrino 
uri and two GeV scale right-handed neutrinos i'r2,3- The keV scale uri has a lifetime 
longer than the age of the universe and is a DM candidate. One of the virtue of this 
model is that the DM particle is already introduced in seesaw mechanism for explaining 
the tiny masses of active neutrinos. 

One interesting implication of uri DM is that uri can interact with Standard 
Model(SM) particles through small mixing with active neutrinos and it might be detected 
in laboratory experiments. However it was shown by some authors |2| that astrophysical 
constraints and the constraints from the mass and mixing of light active neutrinos are very 



strong. Mixing of z/^i with active neutrinos is very small and detecting z/ri in laboratory 
experiment seems difficult. In this article we investigate the mixing of ujn and its density 
in the thermal history of the early universe. We show that relatively large mixing of ujn 
with active neutrinos is allowed in the i/SM when two GeV scale i'r2,3 are degenerate. 
We study the direct detection of the cosmic background of this keV scale u^n DM in beta 
decay experiment. In the following we briefly describe the z/SM and discuss issues related 
to DM. We study the mixing of light neutrinos with heavy neutrinos in the case that two 
GeV scale right-handed neutrinos are degenerate. We analyze the capture of the cosmic 
background of Ufji DM in target of radioactive nuclei ^H and ^''^Ru. 

keV scale z/^^i DM in z/SM: In seesaw mechanism mass matrix of active neutrinos 
is given by 



m,, = V 



Y*M^^Y\ (1) 



'R 



where Mr is the Majorana mass matrix of right-handed neutrinos, Y the matrix of Yukawa 
couplings, V the vacuum expectation value of doublet Higgs in the SM: < H >= (0, v)'^. 
In seesaw mechanism Y can be parameterized as |3|] 

Yv = U{m,f'^P^/^0{M}^Y'\ (2) 

where fhj = diagjm]^ ,"^2 5 ""^3 } ^^^ P^^'^ — diagje*"'^^/^, e*"'^^/^, e*"'^^/^}. Real numbers 
mi{i = 1,2,3) are the masses of light neutrinos and 0j(i = 1,2,3) are Majorana phases. 
(M^)^/^ = diag{(Mj*)^/^, {M^Y^^, (M|)^/^}. U is the neutrino mixing matrix observed in 
neutrino oscillation experiments. O is a complex orthogonal matrix: 0^0 = 00^ = 1. 
Active neutrinos mix with heavy Majorana neutrinos through a mixing matrix 

R = Yv{M*^r\ (3) 

Using Eq. ^ R\s rewritten as 

R = U{m,f'^P^'^0{M*ji)-^'\ (4) 

For convenience we will work in the base that Mi^2,3 are real. 

keV scale vpti can decay to light active neutrinos or photon and its main decay channel 
is the decay to three active neutrinos. For 91 = ^^ l-RnP ~ 10^*^ the lifetime of z/^i 



R 



is estimated ~ 10 s [1] which is much larger than the age of the universe ~ 10 s 
and long enough for a DM candidate, i^m DM can be produced in the early universe 
through active-sterile neutrino oscillation [l|, |j-l6|. For 9f ~ 10^^ significant amount of 
keV scale u^i can be produced. With larger mixing um DM can be over-produced using 
this mechanism. Many other aspects of this model of DM, e.g. constraints from X-ray 
observation and Lyman-g forest, detection of this DM, possible symmetries etc., have 
been analyzed |7l-ll5|. Suggestions beyond the minimal model have also been made and 
other DM production mechanism has been considered 16|-|20|. 

ujii density can be diluted in the early universe and Of is allowed to be much larger if 
significant entropy release is produced in a reheating at MeV temperature scale before the 
Big Bang Nucleosynthesis(BBN). This can be achieved by the decay of a non-relativistic 
particle in the early universe. It's interesting that GeV scale right-handed neutrino is a 
candidate of this non-relativistic particle. However it was noticed that this is a further 
constraint on the seesaw mechanism and large enough entropy production can not be 
obtained because of the constraint from mass and mixing of light neutrinos |2| . In the 
following we show that the analysis in [2] does not really apply to the case when two GeV 
scale right-handed neutrinos are degenerate (M3 = M2) and Of ^ 10~^ is indeed allowed. 

Entropy production by GeV scale ur and neutrino mixing: Using an explicit 
example we show that mass and mixing of light neutrinos and enough entropy production 
by the decay of a GeV right-handed neutrino can all be accommodated in uSM when two 
GeV scale right-handed neutrinos are degenerate. 

Two degenerate right-handed neutrino states can be rewritten in other base using a 
unitary transformation: 



(5) 



where |cp + |sp = 1. The symmetric mass matrix M'^ in this base is no longer diagonal. 
We note that two states |z/^2 3 > have the same mass of |z/r2,3 >• This is just another way 
to write the mass term. For example, if c = 1/a/2 and s = i/\/2, M'^ is obtained (up to 
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a phase factor i ) as 
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(6) 



In Eq. i^ M2 is written in the form of Dirac type. Another way to see this point is to 
note that MrM\^ and M\Mr are not changed by Eq. ^ when M3 = M2. Hence Eq. ^ 
does no change the energy dispersion although Mr is transformed by it. 

When z//j2,3 are degenerate we should transform them to interaction base to understand 
their interaction with active neutrinos. In the interaction base of v'^2 3 ^' i^ obtained from 
Y: Y' = YV'^ where V is the unitary matrix in u'^ = Vvr as shown in Eq. ^. In the 
interaction base Y "^Y' is obtained from Y'^Y to the following form 



' yl vu yh^ 



Y'\Y' 



vtl vi 

,2* n „,2 



(7) 



\2/i3* ^ vij 

Coupling of z/^2 3 with active neutrinos are 2/2,3 times active neutrino mixing. 

As an example, we consider normal mass hierarchy and matrix O of the following form 



O 



1 

ce S0 
-S0 ce 



where ce = cos 6, se = sin 6 and 9 = x + iy. x and y are real numbers. Diagonalizing the 
second and third entries of Y'^Y using Eq. ([5]) we obtain c = cos a, s = sin a e*^ where 

2D 



tan 2a 



A-B'' 



De"" = {m24se - m^slce)Ml'^Ml'\ 



rl/2»^l/2| 



and D = \{m2C*gSg - m^slce)M^''' M^'''\, A = M2{m2\ce\^ + m^lsel^), B 
mslcep). Three eigenvalues are found: yl = rriiMi/v'^, 



yl 



ICIv' 



A + 5 + v/(A + 5)2-4C 



1/3 






(9) 
(10) 

M3(m2|se|2 + 

(11) 
(12) 



where C = m2m^M2M^. 

For \y\ > 1 it's easy to see that \ce\ ^ Is^l ^ |el^l, A ^ \M2{m2 + m3)e^l^', 5 ^ 
jM^{m2 + m3)e^l^l and we find 

^' (m2 + m3)(M2 + M3) ' ^ ^ 

yl ~ i^(^2 + m3)(M2 + M3)e2|^l (14) 

For normal mass hierarchy we have m^ "^ yAmf^ ^ 0.05 eV, m2 ~ \/'^^oiar ~ 0.009 
eV and mi ^ m2. Using M3 = M2 we get 



2m2M2 1^1 m3M2 21,1 



^,^^^^^^_,,, ^3^^^^^,,_ ^^^^ 



In Eq. ( TT5I) we see that couphng of z/^2 with active neutrinos are suppressed by e^'^' if 
\y\ is large and coupling of z/^3 is enhanced. We note that cos 2a; — > 0, \f3\ — > 7r/2 in the 
limit li/l — )■ 00 and the mass matrix Mjj approaches the mass matrix shown in Eq. ([6]). 
The mixing R' in the interaction base can be similarly obtained. We obtain 6^2 3 = 



12,13\ ^S 



e[' = ^, (16) 

^ Ml' ^ ' 
^, 2m^^_2^y^ 

2 M2 ^ ^ 

0;2^J!^e2H. (18) 



^ 2M9 



The decay rate of z/^2 (^^ ^2 < ^w) is [2'] 



r = ^IMi.^'2 ^ g^^2M| ^^.,|,| 



1927r3^ "" ' 967r3 ^ ^ ' ^^^^ 

where a factor 2 has been included to account for the charge conjugation processes and 
a factor 1/2 has been included to account for the Majorana nature of i^^2- -^ is a factor 



a 



which accounts for effects of various final states J2|. F ?^ 16.0 for rrih ^ M2 < M^. The 
lifetime of z/^2 is 

e^ls^l 0.01 eV /30 GeV\^ , , 

T ^ 0.44 s X . 20 

(2000)2 7712 \ M2 J 

It is worth pointing out that u'^2 can not decay through oscillation to z/^3 because the 
probability of oscillation to z/^3 vanishes when z/^2 3 ^'^s degenerate. The above analysis 



and the later discussions can also be applied to quasi-degenerate i'r2^2, if their mass dif- 
ference is so small that no significant oscillation can happen before v'^2 decays. For our 
whole discussion to be valid for quasi-degenerate i'r2,^^ it's sufficient to assume that no 
significant oscillation can happen before the BBN time which is ~ 1 s. This condition is 
IM2 - M3I X 1 s < n, or IM2 - M3I « 6.58 x lO'^^ j^eV. 

Assuming that z/^2 is in thermal equilibrium at very high temperature and decouples 



when it is relativistic, the entropy release produced by the decay of z/^2 i^ 1^1 



S ^ 0.76- 



2 g^VT'Mfi 
'0.01eV\^/V30 GeV\ e'^l 



_ /0.01eV\'^V30 GeV\ 

^ 96 X X , (21) 

\ m2 J \ M2 J 2000' ^ ^ 

where S = Sf/si. Sf and Sj are entropies just after and before the decay of 1/^2- fi'Af = 2 is 

the number of degrees of freedom of z/^2 ■ 9* is the effective number of degrees of freedom 

at z/^2 freeze-out which is assumed at electro- weak scale, g^ ~ 101.5 when excluding top 

quark, g^ ^ 10.75 is the effective number of degrees just after the entropy release. 

In Eq. ( I2T]) we can see that S ~ 100 can be achieved for e'^' ~ 2000. According to 

Eq. ( I20I) . for this range of parameter space z/^2 ^^s lifetime ^ 1 s. So it decays before 



the BBN and does not spoil the prediction of the BBN [22|. It's also easy to see that 
without e'^' factor M2 ^ 1 GeV is needed to make S larger than 100. However, M2 <^ 1 
GeV is not allowed because it would give a too long lifetime to z/^2 ^^^ i^ would spoil the 
prediction of the BBN. We emphasize that previous negative conclusion on the entropy 
production by iyR2 [21] does not apply to the case we are considering. They worked in 
the base that mass matrix is diagonalized and did not pay attention to the fact that we 
should analyze the interaction of GeV scale right-handed neutrinos in their interaction 
base when they are degenerate. 

We point out that in our example 6^ ~ 10~^ can be obtained by taking mi ~ 10^'^ 
eV which is allowed by the neutrino oscillation experiments. Astrophysical observation of 



the decay uri — )■ z/7 is able to constrain 6^ 231] 



.?<i,8x:o-»(iM)^. 



(22) 



We see that Mi ^ 2 keV is needed for 6-^ reaching ~ 10 ®. Mi is constrained by the 
observation of dwarf spherical galaxies: Mi ^1—2 keV pLL] • Mi is also constrained by 



Lyman-a forest. In the case we are considering the Lynian-a bound is Mi ^1.6 keV when 
S ~ 100 [2] which is re-scaled by a factor S^^^^ compared to the bound for non-resonant 



O 



(23) 



production of ujn DM [lOj. We further note that X-ray and Lyman-a 10| constraints 
become weaker when upa DM accounts for part of the DM in the universe, e.g. ~ 40% of 
total DM energy density. In summary Mi ^ 2 keV is allowed by the present constraints 
and 9^ is allowed to reach ~ 10~^ |2|. 

We note that z/^2 has very small Yukawa coupling and can not come into thermal equi- 
librium due to this interaction. Other interaction of z/^2 is needed to make it populated at 
temperature of electro-weak scale. Detailed model for it will be explored in future works. 
The example shown is a particular case with normal hierarchy. The quasi-degenerate 
mass pattern of light neutrinos would give too large mixing of keV scale right-handed 
neutrino with active neutrinos and is not compatible with the i/SM [l|, |l5| . For inverted 
mass hierarchy, we can choose matrix O as 

' cos 6 sin 9 * 

— sin 9 cos 9 
yl y 

Similar results can be obtained in the case that two GeV scale i'r2,3 are degenerate. 

u^i capture by radioactive nuclei: u^i can interact with SM particles through 
mixing with active neutrinos. The four-Fermion interaction of ujn with electron is 

C 

^L = -^RlinYiQA - 759 a)p 

xz>Ri7^(l-75)e + /i.c., (24) 

where n and p stand for neutron and proton. gv,A are vector and axial- vector form factors. 
According to Eq. (IMl) uri DM in the universe can be captured by radioactive nuclei in 
processes 

uri + A ^ B~ + e~^ or um + A ^ B'^ + e~ . (25) 

The energy of electron produced in this process is 

Ee = me + Qp + Ml, (26) 

where Q/j is the end point kinetic energy of the beta decay {A -^ B^ + e^ + Ue) or anti-beta 
decay {A — t- B~ + e"^ + z/g). For Mi <^ Q^, the cross section of this process is similar to 



8 

that of the capture of the cosmic rehc neutrinos 2J, |25|. The capture of um DM produces 
mono-energetic electrons beyond the end point of the spectrum of beta decay or anti-beta 
decay. For super-allowed transition the cross section is 

O-captVu = -r. PeEeF{Z, Ee)\Rel\'^, (27) 

7*1/2 



where pe = \/2me{Qji + M\) and -F(2', ii^g) is the Fermi function. In Eq. fl27|) a factor 
of 1/2 arising from the Majorana nature of vji\ has been included. The cross section is 
normalized to the beta decay rate using ti/2, the half-life of A nucleus and the factor / 



/=/ F{Z,E'MKKPudE',- (2^ 

J m,. 



Pg, E'^ are momentum and energy of electron (or positron ). p'^and E^ are the momentum 
and energy of anti-neutrino (or neutrino) in beta decay (or anti-beta decay). 

Consider uri capture by Tritium: uri + ^H -^ ^He + e^ which has Q/s = 18.59 keV 
and ti/2 = 12.3 year. Using parameters in [24] we find that for Mi <^ Qp the cross section 
of this process is (XcaptVu ~ 3.9 x 10^^^ cm^^ x c x |-ReiP where c is the speed of light. The 
event rate, (JcaptVunujii, is estimated 

N ^ 0.71 year~i x ^'-m J^eil^ 'H 

^ 105 cm~3 10-6 10 kg ^ ^ 

where ujn DM is assumed to account for a significant part of DM and its number density 

is estimated 

In Eq. f pOj) we have used reference density 0.3 GeV cm^^, the estimated DM density 
in the galactic halo at the position of the solar system. This estimation of local density 
is extrapolated from the astrophysical observation of Milky Way and the model of DM 
halo which does not depend on whether the DM is warm or is cold. This estimation is 
applicable to warm DM. 

We find that ^"^Ru, which has Qp = 39.4 keV and ti/2 = 373.6 days, is also good to 
detect Ufii. The cross section of uri capture by ^"^Ru is (JcaptVu ~ 2.94 x lO'^^^ cm"^ x 
c X l-ReiP for Ml <^ Qj3. The event rate of uri capture is 

A^ ^ 16 year-^ x /^^ ' ^^' . 31 

■^ 105 cm-3 10-6 10 Ton ^ ^ 



We note that the production rate of ^H in reactor is 0.01% and the production rate of 
^°^Ru is 0.4%. If 10 kg Tritium are produced per year in reactors, around 12 Tons of 
^°^Ru are produced. We note that the reference number Eq. (!30|) used in Eq. (!3T]) is a 
conservative estimate of the local number density. It's possible that the solar system is 
located in a sub-halo in which the local DM density is several orders of magnitude larger 
than the galactic value. If this happens, capture rate is several orders of magnitude larger. 
We emphasize that the events of ujn capture are mono-energetic electrons which have 
energy well separated from the (3 decay spectrum. The background of the (3 decay events 



does not affect the detection of ujn. Moreover, 



;his experiment does not require mea- 



surement of very high precision as in KATRIN |27| | or in detecting cosmic background 



neutrinos |24J . |25| . Rather than using gaseous Tritium in KATRIN, large volume of ^H 
or ^°^Ru target in solid state can be used in ujn capture experiment. This experiment 
might be done in the near future. We note that future X-ray observation may improve 
the constraint on 6^ which will reduce the expected event rate of z/^i capture. However 
the scenario of keV scale right-handed neutrino as DM candidate is hard to be ruled out 
and the DM capture by radioactive nuclei proposed in this article is one way to detect 
this keV scale DM candidate. 

We note that background events of the capture of solar pp neutrinos are negligible. 
According to the standard solar model the flux of pp neutrinos of energy ^ 10 keV is 
about 8. X 10^ cm~^ s"^ |26| and the corresponding density of these pp neutrinos is about 
2.7 X 10"'* cm~^. The event rate of the capture of pp neutrinos of energy ^ 10 keV is 
~ 4.0 X 10"^ year"^ for 10 kg of Tritium and is ~ 8.5 x 10"^ year^^ for 10 ton of ^°^Ru. 
Effects of low energy solar neutrinos can be neglected in discussing the capture of um 
DM. 

Recent analysis of X-ray observations of local dwarf Willian 1 show evidence that the 



Z//J1 DM may have mass around 5 keV with mixing |-RnP ~ 10^^ 28|]. Another analysis 
of the X-ray observation of the galactic center suggests that z/^i DM has mass around 17 
keV with mixing |-R;iP ~ 10~^^ 29|]. This small neutrino mixing would give too small 



event rates of um capture by radioactive nuclei unless we are living in a sub-halo of DM. 
Conclusion: In conclusion we have considered several issues of the keV scale right- 
handed Z//J1 DM in the z/SM. We have shown that a GeV scale right-handed neutrino in 
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the z/SM can have sufficient small couplings to active neutrinos and its decay can lead 
to large amount of entropy release ( ~ 100 ) at MeV temperature scale. This happens 
when two GeV scale right-handed neutrinos are degenerate. Two degenerate right-handed 
neutrinos in their interaction base can have very different strengths of couplings with 
active neutrinos and a large suppression to the Yukawa couplings of one of the GeV scale 
vrS can be achieved. Mass and mixing of light neutrinos can be correctly explained in 
the model considered. Density of vri DM can be diluted by a factor larger than 100 in the 
reheating produced by the decay of this GeV scale right-handed neutrino. The mixing of 
vri with active neutrinos is allowed to be as large as 9\ ~ 10 ~® for Mi k. 2 keV, reaching 
the bound from X-ray observation. 

We have discussed the capture of the cosmic background of keV scale vri DM by 
radioactive nuclei ^H and ^*^^Ru. vri capture produces events of mono-energetic electrons 
with keV scale energy away from the beta decay spectrum. The background of beta decay 
events does not affect the detection of vri capture events and this experiment does not 
require measurement of very high precision as in experiment of direct measurement of 
neutrino mass or in detecting cosmic background neutrinos. Rather than using gaseous 
Tritium in experiment of direct measurement of neutrino mass, large volume of "^H or 
^°^Ru target in solid state can be used in vm capture experiment. We find that there are 
about 0.7 events per year on 10 kg tritium target and about 16 events per year on 10 Ton 
^°^Ru target for mixing |-ReiP = 10^^ and Uj,^^ = 10^ cm^^. 

We comment that if we are living in a sub-halo of DM the event rate can be much 
larger. Two examples with ^H and ^°^Ru targets are related to beta decay. It is 
interesting if suitable radioactive nucleus can be found to do vri capture in anti-beta 
decay experiment. We emphasize that detecting vri DM on radioactive nuclei is a 
new type of experiment for detecting DM in the universe. It is very interesting if 
other radioactive nucleus, better than ^H and ^'^^Ru, can be found. This radioactive 
nucleus should have reasonable lifetime, large enough capture rate and reasonably large 
production rate in reactors or is available in nature. It's worth further exploration. 
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